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Abstract 

We review recent advances in modeling global-scale convection and 
dynamo processes with the Anelastic Spherical Harmonic (ASH) code. 
In particular, we have recently achieved the first global-scale solar con- 
vection simulations that exhibit turbulent pumping of magnetic flux 
into a simulated tachocline and the subsequent organization and am- 
plification of toroidal field structures by rotational shear. The presence 
of a tachocline not only promotes the generation of mean toroidal flux, 
but it also enhances and stabilizes the mean poloidal field throughout 
the convection zone, promoting dipolar structure with less frequent po- 
larity reversals. The magnetic field generated by a convective dynamo 
with a tachocline and overshoot region is also more helical overall, with 
a sign reversal in the northern and southern hemispheres. Toroidal 
tachocline fields exhibit little indication of magnetic buoyancy insta- 
bilities but may be undergoing magneto-shear instabilities. 



1 Introduction 

There is little doubt that turbulent convection in the solar envelope must 
convert kinetic energy into magnetic energy by means of hydromagnetic 
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dynamo action. However, theoretical arguements and numerical simulations 
of turbulent dynamos suggest that much of this field must be in a disordered 
state, dominated by small-scale stochastic fluctuations. The challenge of 
global solar dynamo theory is to account for the ordered component of the 
magnetic field, that which gives rise to sunspots, large-scale fields features 
such as dipolar and quadrupolar moments, and the 22-year activity cycle. 

Advances in solar dynamo theory over the past several decades suggest 
that boundary layers may play a crucial role in generating this ordered mag- 
netic field component. In particular, the solar tachocline located near the 
base of the convective envelope may be a breeding ground for the manufac- 
ture and accumulation of large-scale magnetic flux. There are several reasons 
for this. (I) Convective expulsion, turbulent pumping, and global meridional 
circulations continually transport magnetic field into the tachocline region 
where it can accumulate. (II) Rotational shear amplifies mean toroidal fields 
through the O-effect and suppresses smaller-scale fields by enhancing ohmic 
diffusion, thus promoting the generation of strong, ordered toroidal flux. 
(Ill) The subadiabatic stratification inhibits magnetic buoyancy instabil- 
ities, enabling long-term storage of magnetic flux. (IV) The presence of 
strong, stable large-scale toroidal flux and electrical current structures in 
the tachocline may provide a magnetic inertia, enhancing and stabilizing 
the poloidal held component as well. Other processes may also contribute 
to the generation of ordered magnetic fields and cyclic activity, including 
the breakup of active regions in the upper boundary layer (the Babcock- 
Leighton mechanism) and advection of magnetic fields by global meridional 
circulations (e.g. Charbonneau 2005). 

Here we review insight into the nature of global solar and stellar dynamo 
processes gained from 3D MHD simulations of turbulent convection in ro- 
tating spherical shells. We focus on two representative dynamo simulations 
in a solar context, one in which the computational domain is limited to the 
convective envelope and one that incorporates convective overshoot into an 
underlying radiative zone where uniform rotation is imposed, thus creating a 
mock tachocline. In §2-3 we compare the two simulations and assess how the 
presence of this tachocline and overshoot region alters the global dynamo. 
The most apparent difference is the presence of strong, organized toroidal 
fields in the tachocline region of the latter simulation. In §4 we discuss how 
such fields may be maintained and address their stability and evolution. In 
§5 we briefiy summarize the implications of these and other stellar dynamo 
simulations. 
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Figure 1: (a-c) Illustration of convective patterns and mean flows in Case 
B. (a) Orthographic projection of the radial velocity, Vr near the middle 
of the convection zone (r = 0.84i2), with the north pole titled 35° toward 
the line of sight. Dashed lines indicate the equator and several meridians. 
Bright tones (yellow/orange in online version) denote upflow and dark tones 
(blue/black) denote downflow. (6) Angular velocity and (c) meridional cir- 
culation averaged over longitude and time (90 days), the latter represented 
as streamlines of the mass flux. The dashed line indicates the base of the 
convection zone, (d) The angular velocity in Case A is also shown for com- 
parison. Bright (pink/red) and dark (blue/black) tones denote faster and 
slower rotation in (b) and (d), ranging from 350 to 450 nHz. Bright (red) and 
dark (blue) tones in (c) denote clockwise and counter-clockwise circulations 
respectively, with the flow speed ranging from 2-20 m s~^. 



2 Simulations and Flow Features 

In order to assess the role of a tachocline in global dynamo action we consider 
two 3D MHD simulations of global solar convection. Both are sustained 
dynamos that extend over multiple ohmic diffusion time scales. The first is 
described in detail by Brun, Miesch &; Toomre (2004; hereafter BMT04), and 
is referred to there as Case M3. Here we refer to it as Case A. The second, 
which we refer to here as Case B, is described further by Browning et al. 
(2006; hereafter BMBT06; see also Browning et al. 2007). Both simulations 
are based on the Anelastic Spherical Harmonic (ASH) code with solar values 
used for the luminosity, rotation rate, and background stratification. 

The principle difference between Cases A and B is the presence of a 
tachocline in the latter. Whereas the computational domain in Case A 



4 



Miesch et al. 



stops at the base of the convection zone, r = 0.72R where R is the solar 
radius, Case B incorporates a portion of the subadiabatic radiative interior, 
extending downward to r = 0.62i?. Both simulations stop short of the 
solar photosphere, with upper boundaries at (A) 0.97-R and (B) 0.96R, in 
order to avoid complications in the photospheric boundary layer such as 
granulation, ionization, and the transition to radiative energy transfer which 
cannot be easily accounted for in our global anelastic modeling framework. 
The kinematic viscosity, v, is somewhat higher in Case B, 2 x 10^^ cm^ s^^ 
in the mid convection zone versus 1.4 x 10^^ cm^ in Case A, but the 
magnetic diffusivity rj is somewhat lower {Pm = 4 in Case A and = 8 
in Case B where Pm = '^/v)- Both simulations have comparable spatial 
resolution. For further details see BMT04 and BMBT06. 

A tachocline is imposed in Case B by means of a drag force in the ra- 
diative interior that imposes uniform rotation. This is intended to take into 
account tachocline confinement processes that are not captured by the sim- 
ulation due to limited spatial resolution and temporal coverage. A thermal 
forcing is also applied near the base of the convection zone, warming the 
poles relative to the equator. This is intended to help sustain a strong rota- 
tional shear across the tachocline by means of baroclinic torques that par- 
tially offset the sub-grid-scale diffusion and that compensate for the extent 
of the overshoot region which is wider than in the Sun due to computational 
limitations. We emphasize that the rotational shear in the convective en- 
velope is still maintained by the convective Reynolds stress; the mechanical 
and thermal forcing merely promotes a sharp tachocline. Such an approach 
cannot provide much insight into tachocline confinement processes but it 
can be used to investigate the influence of a tachocline on global dynamo 
action, which is our motivation. For further discussion see BMBT06. 

The convective structure in Cases A and B is similar, dominated at low 
latitudes in the mid convection zone by columnar convection cells aligned 
with the rotation axis (Fig. [Hi) , as in comparable non- magnetic simulations 
(reviewed by Miesch & Toomre 2009). Near the surface both exhibit an 
interconnected network of downflow lanes laced by vertical magnetic flux 
(BMT04). 

The differential rotation profile in both cases is roughly solar-like (Fig.[TI), 
d), with a monotonic decrease in angular velocity with latitude. However, 
both exhibit more cylindrical alignment than the solar rotation profile as 
revealed by helioseismology. Furthermore, the angular velocity contrast AO 
between equator and pole is much weaker in Case B (about 13% compared 
to 36% in Case A). This is typical of penetrative convection simulations and 
may arise in part from our artificially wide overshoot region (Miesch 2007a) 
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Figure 2: Magentic (black; red in online version) and kinetic (grey; blue) 
energy density in (a) Case A and (b) Case B, averaged over horizontal sur- 
faces and over time (4 days). Solid lines denote the energy in the mean 
toroidal field and the differential rotation, dotted lines that in the mean 
poloidal field and the meridional circulation, and dashed lines that in the 
non-axisymmetric field and flow components. The vertical dotted line in (6) 
denotes the base of the convection zone, (c) Molleweide projection of the 
longitudinal field component in the tachocline of Case B (r = 0.67i?) at 
a time when m = 1 structure is particularly evident. Bright and dark tones 
denote eastward and westward fields respectively. 



and in part from Lorentz forces associated with stronger mean fields (§3). 
Despite the weaker AQ, Case B possesses a tachocline; substantial radial 
shear near the base of the convection zone maintained by our mechanical 
and thermal forcing. 

The most prominent feature of the mean meridional circulation in Case B 
is an equatorward flow at the base of the convection zone with an amplitude 
of about 8 m s~^ (Fig. This may in part be attributed to the thermal 
forcing but other simulations without such forcing also exhibit equatorward 
circulation near the base of the convection zone, albeit somewhat weaker 
(several m s~^; see Miesch et al. 2000). This has important implications for 
flux-transport dynamo models as we address in §4. The circulation proflle 
in the bulk of the convection zone is multi-celled, as in Case A (BMT04). 



3 Mean and Fluctuating Fields 

As discussed in §1, turbulent flows beget turbulent flelds. In both Cases A 
and B, dynamo action produces strong fluctuating (non-axisymmetric) flelds 
that account for 95-98% of the total magnetic energy in the convection zone 
(Fig. [2]3, b). Both mean and fluctuating fleld components are somewhat 



6 



Miesch et al. 




Figure 3: Mean (a, b) poloidal magnetic fields and (c, d) toroidal magnetic 
fields in (a, c) Case A and (6, d) Case B, averaged over longitude and time 
(90 days). White (red in online version) tones denote a clockwise poloidal 
field orientation and eastward toroidal field whereas black (blue) tones de- 
note the opposite. Saturation levels for the poloidal and toroidal color ta- 
bles are ±20G and ib2000G respectively. Poloidal fields in (a, b) include a 
potential-field extrapolation out to 1.5-R and the base of the convection zone 
is indicated by a dashed line in (6, d). 

stronger in Case B, due in part to the lower magnetic diffusivity (§2). More 
remarkably, the magnetic energy in the mean toroidal field of Case B rises 
steadily toward the base of the convection zone, peaking in the tachocline 
where it reaches equipartition with the fluctuating field components and also 
with the fluctuating kinetic energy. Note that the drop in the differential 
rotation kinetic energy below r = 0.69i? is due primarily to the mechanical 
forcing (§2). 

The structure of the tachocline field in Case B is dominated by strong 
(2-5 kC) toroidal bands or extended magnetic layers that are roughly anti- 
symmetric about the equator. These toroidal bands are generally axisym- 
metric as illustrated in Figure [3li but they occasionally exhibit a pronounced 
m = 1 component as demonstrated in Figure [2]:;, where m is the longitudinal 
wavenumber. This m = 1 behavior may be indicative of magneto-shear in- 
stabilities as we discuss in §4. The toroidal bands are remarkably persistent, 
remaining intact for thousands of days, punctuated sporadically by brief in- 
tervals (< 200 days) of roughly symmetric parity about the equator after 
which the antisymmetric parity is re-established with the same polarity as 
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Figure 4: Time evolution of the dipole (black) and quadrupole (grey; red in 
online version) moments in (a, c) Case A and (6, d) Case B. The top row 
corresponds to the outer boundary in each case (a) 0.97i? and (b) 0.96i? 
whereas the bottom row corresponds to the base of the convection zone in 
Case A (c) 0.72R and to the underlying tachocline in Case B (d) 0.67R. 
Points denote instantaneous snapshots and lines denote running six-month 
averages. 



before (BMBT06, Browning et al. 2007). The simulation has spanned over 
two decades and has not yet undergone a sustained polarity reversal of the 
mean toroidal field in the tachocline. 

In contrast to the strong, persistent bands in the tachocline, the mean 
toroidal field in the convection zone of both cases A and B is relatively weak 
(~ 1 kC), disordered, and transient (Fig. [S]:, d). This lends support to the 
arguments put forth in §1 on why the tachocline might be a prime place to 
generate large-scale magnetic flux. The mean poloidal field in Case B (Fig. 
[3^) also appears stronger and more ordered than in Case A (Fig. [3h), with 
a more prominent dipole moment. This is demonstrated further in Figure H] 
which shows the dipole and quadropole moment in each case over thousands 
of days of evolution. 
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Figure 5: The fluctuating current helicity J'-B' is shown for (a) Case A 
and (6) Case B averaged over longitude and over the same 90-day time 
interval as in Figures [1] and [3l White and black (red and blue in online 
version) denote positive and negative values and the color table saturates 
at ±2 X lO''' statamp cm~^ G. (c) Kinetic helicity in Case A (black line) 
and Case B (grey; blue), averaged over the northern hemisphere and time 
(4 days). Dashed lines in (6) and (c) denote the base of the convection zone 
in Case B. 

In Case A, the dipole and quadrupole moments are comparable in am- 
plitude and erratically reverse sign on a time scale of roughly 500 days. By 
contrast, the dipole moment in Case B is stronger both in absolute mag- 
nitude and relative to the quadrupole component, although this is less so 
toward the end of the interval shown in Figure [H The dipole moment is also 
more stable in Case B, maintaining the same sign through most of the simu- 
lation apart from a brief reversal near the top of the shell at t ~ 6000 days. 
As with any numerical simulation, however, these conclusions are based on 
a limited time interval. Figure [Hi in particular suggests that the dipole 
moment may be waning over the course of several decades. Further time 
evolution is necessary to clarify any potential long-term trends, including 
the possiblity of systematic polarity reversals. 

Another remarkable feature of the magnetic field in Case B is illustrated 
in Figure [5f>; it is helical. Shown is the current helicity in the fluctuating 
field component, H'^ = J'-B', where J = (c/47r)VxS is the current density 
and primes indicate that the longitudinal mean has been subtracted out. 
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Although the magnetic hehcity Hm = A-B is of greater theoretical interest 
(where A is the vector potential, defined such that B = VxA), He has 
the advantage that it is guauge invariant and is more readily observable. 
Both He and Hm quantify the degree to which the magnetic field topology 
is twisted and is not reflectionally symmetric. 

Magnetic helicity has profound implications for large-scale dynamo ac- 
tion. The turbulent a-effect, whereby large-scale fields are generated by 
small-scale velocity fields, is intimately tied to the upscale spectral transfer 
of magnetic helicity. If the transfer is local, this is referred to as an inverse 
cascade. However, the ensuing buildup of small-scale helicity can suppress 
the inverse cascade if it is not dissipated or removed through the bound- 
aries. This suppression can be so dramatic that it has been referred to as 
catastrophic a-quenching and if it occurs, it would imply that the turbulent 
a-effect cannot account for the mean fields observed in the Sun and other 
stars. For further discussion see Brandenburg & Subramanian (2005) and 
Miesch & Toomre (2009). 

The helical nature of the fiuctuating magnetic field in Case B (Fig. [5f)) 
suggests that the presence of a tachocline and overshoot region may miti- 
gate catastrophic a quenching. The convection may be shunting small-scale 
helicity into the stable zone where it is stored or dissipated, thus promoting 
helicity generation, and possibly self-organization processes associated with 
upscale helicity transfer, in the convection zone. By contrast, the fluctuat- 
ing field in Case A (Fig. [5jz) is not particularly helical and does not exhibit 
the asymmetry about the equatorial plane evident in Case B. 

The sign of the current helicity in Case B is positive in the north and 
negative in the south, which is opposite to that inferred from photospheric 
and coronal measurements (e.g. Pevtsov & Balasubramaniam 2003). How- 
ever, such measurements are concerned with field structures such as active 
regions that are associated with the large-scale field component. If such 
fields are generated through an upscale transfer of magnetic helicity, one 
would expect that the sense of the helicity should be opposite to that in the 
fluctuating component. Indeed, the negative apparent in the northern 
hemisphere in Fig. [3li coupled with the positive dipole moment apparent 
in Fig. [3li implies a negative current helicity for the mean fleld component. 
Thus, the sign of the current helicity in Case B may be consistent with solar 
observations. 

There are many caveats to this simplifled picture, including the poorly 
understood details of how helicial flux tubes form in the tachocline and 
emerge. Nevertheless, it is clear that the dynamical aspects of magnetic he- 
licity offer a fertile ground for further insight into the subtleties of global so- 



10 



Miesch et al. 



lar and stellar dynamos, particularly with regard to influence of a tachocline. 

In contrast to H'^, the kinetic helicity, = u-v where v and u are the 
fluid velocity and vorticity, is similar in Cases A and B (Fig. [5];) . The Hk 
profile in both cases is antisymmetric about the equator, with a negative 
sign in the northern hemisphere through most of the convection zone. Near 
the base of the convection zone the sign of Hk reverses due largely to the in- 
fluence of the Coriolis force on diverging downflows. Kinetic helicity profiles 
such as this are a well known feature of rotating compressible convection 
(e.g. Oilman 1983; Miesch et al. 2000). 

4 Maintenance and Evolution of Tachocline Fields 

According to our current understanding of the solar dynamo, we may expect 
that toroidal fields generated in the tachocline by the Q-efFect should mi- 
grate equatorward, intermittently spawning flux tubes as they go which then 
emerge through the photosphere as active regions. In this section we briefly 
address how our simulations relate to this prevailing dynamo paradigm. 

To this end, we must first understand how the toroidal fields in the 
tachocline of Case B are established and maintained. A detailed analysis 
will await a future paper but here we can confirm that the principle source 
term in generating mean toroidal flux in the tachocline is indeed the r^-effect 
arising from radial angular velocity shear dQ/dr, as might be expected based 
on the prevailing paradigm. The latitudinal shear dfl/dO also contributes, 
as does the radial component of the turbulent emf: E'^ = v'qB'^ — v'^B'g. 

These generation terms work to offset radial diffusion and the latitudinal 
advection of toroidal flux by the meridional circulation which both act to 
erode the tachocline field (latitudinal diffusion is negligible) . The latitudinal 
component of the turbulent emf, Eg = v'^Bl — v'^B'^. also acts as a sink term 
on average, at least over the 6930-day time interval shown in Fig. Uf), d. 
This is somewhat surprising since turbulent pumping of toroidal flux, if it 
occurs, would be reflected in the E'q term. Thus, turbulent pumping acts to 
erode the tachocline fleld in Case B rather than maintain it. Although this 
may be contrary to theoretical expectations, it is still consistent with the 
concept of turbulent pumping since the mean toroidal flux in the convection 
zone of Case B is opposite in sign (albeit weaker in magnitude) relative to 
the tachocline field. 

Conventional wisdom suggests that as toroidal flux is amplified by rota- 
tional shear in the tachocline it will eventually become susceptible to mag- 
netic buoyancy instabilies, forming isolated flux tubes which then rise toward 
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the solar surface. Although some strong, transient toroidal ribbons in Case 
B do exhibit a slight density deficit which may induce them to rise, the per- 
sisent magnetic layer evident in Figure [3]i appears to be largely stable. The 
proximate reason for this is because the layer does not exhibit any significant 
density evacuation relative to its surroundings but the ultimate explanation 
is still under investigation. It may be that our subgrid-scale diffusion and 
our artificially wide overshoot region limits the radial magnetic field gradi- 
ents that can be achieved. In any case, it is clear that turbulent flows and 
fields in Case B contribute to the mechanical equilibrium of the magnetic 
layer and furthermore, the layer is not an isolated magnetic surface; rather, 
it has fleld line connectivity to the entire convection zone. Thus, conditions 
here are far from the idealized equilibrium states typically considered in 
investigations of magnetic buoyancy instabilities. 

Toroidal magnetic fields in the tachocline are also susceptible to another 
type of instability fed by latitudinal shear, as described by Gilman & Fox 
(1997). Although higher wavenumber modes may also occur, the most vig- 
orous and robust mode for broad toroidal field profiles (antisymmetric about 
the equatorial plane) is the m = 1 clam-shell instability whereby toroidal 
loops tip out of phase, reconnecting across the equator at one meridian and 
spreading poleward at the antipode. If the rotational shear is maintained 
and if the poloidal fleld is continually replenished, such clam-shell instabili- 
ties can occur indeflnitely (Miesch 2007b). 

The presence of sustained clam-shell instabilities in Case B is suggested 
by the prominent m = 1 structure in snapshots such as that shown in Fig. 
[2I3. If such instabilites are indeed occurring, one would expect to see a quasi- 
periodic energy exchange between the m = and the m = 1 components of 
the toroidal fleld (Miesch 2007b). We will address this in a future paper. 

In mean-fleld dynamo models, the equatorward propagation of toroidal 
flux implied by the solar butterfly diagram is usually attributed to one of 
two mechanisms (that are not mutually exclusive). The first is advection 
by an equatorward meridional circulation near the base of the convection 
zone and the second is the propagation of a dynamo wave induced by rota- 
tional shear and turbulent fleld generation (the a-effect). In the latter case, 
equatorward propagation is achieved if the product H/^dCl/dr > in the 
northern hemisphere (e.g. Charbonneau 2005). 

Both of these propagation mechanisms could in principle be occurring in 
Case B. At mid and low latitudes near the base of the convection zone there 
is a strong equatorward meridional circulation (Fig. [T]^) and Hi^ and dVt/dr 
are both positive in the northern hemisphere (Figs. [H), [Sj^). However, the 
amplitude of the circulation and the kinetic helicity peak above the principle 
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toroidal flux concentration. Perhaps not surprisingly, then, there is little 
indication for any latitudinal propagation of the tacholine field in Case B. 

5 Conclusion 

Although 3D MHD dynamo simulations cannot capture all processes of rel- 
evance to the global solar dynamo, they can provide crucial insight into key 
ingredients. In particular, understanding the subtle interaction between tur- 
bulent field generation in the convection zone and the generation of toroidal 
flux by rotational shear in the tachocline stands to benefit greatly from 
high-resolution MHD simulations of penetrative convection. By providing a 
reservoir for expelling, amplifying and storing magnetic flux, the tachocline 
influences the global behavior of the dynamo, promoting strong, stable mean 
fields and helical magnetic topologies throughout the convection zone. 

Throughout this brief paper our emphasis has been on the Sun but stel- 
lar dynamo simulations have also produced substaintial insights in recent 
years. Highlights include the generation of strong toroidal flux structures 
in rapidly-rotating solar-likc stars (Brown et al. 2007), enhancement of core 
dynamo action in A stars by fossil envelope fields (Featherstone et al. 2007), 
quenching of AO by the Lorentz force in the deep convective shells of A 
and M stars (Brun, Browning &; Toomre 2005; Browning 2008), and the 
saturation of Afi and mean field strengths at high 17 (Christcnsen & Aubert 
2006). The latter was done in planetary context but may also have impor- 
tant implications for stars (sec Miesch &: Toomre 2009). 

We thank Nicholas Featherstone, Kyle Augustson and Nicholas Nelson 
for numerous discussions on all aspects of MHD dynamo simulations and 
Keith MacGregor for helpful comments on the manuscript. The work pre- 
sented here was supported by NASA through the Heliophysics Theory Pro- 
gram grant NNG05G124G. The simulations were carried out with NSF PACI 
support of PSC, SDSC, NCSA, NASA support of Project Columbia, and 
through the CEA resource of CCRT and CNRS-IDRIS in Prance. 
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